How do the small planetary satellites rotate? 
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Abstract 

We investigate the problem of the typical rotation states of the small planetary 
satellites from the viewpoint of the dynamical stability of their rotation. We show that 
the majority of the discovered satellites with unknown rotation periods cannot rotate 
synchronously, because no stable synchronous 1:1 spin-orbit state exists for them. They 
rotate either much faster than synchronously (those tidally unevolved) or, what is much 
less probable, chaotically (tidally evolved objects or captured slow rotators). 
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The majority of planetary satellites with known rotation states rotates synchronously 
(like the Moon, facing one side towards a planet), i.e., they move in synchronous spin-orbit 
resonance 1:1. The data of the NASA reference guide (http: / / solarsystem.nasa.gov/planets/) 
combined with additional data (Maris et al., 2001, 2007; Grav et al., 2003) implies that, of 
the 33 satellites with known rotation periods, 25 rotate synchronously. 

For the tidally evolved satellites, this observational fact is theoretically expected. The 
planar rotation (i.e., the rotation with the spin axis orthogonal to the orbital plane) in 
synchronous 1:1 resonance with the orbital motion is the most likely final mode of the long- 
term tidal evolution of the rotational motion of planetary satellites (Goldreich and Peale, 
1966; Peale, 1977). In this final mode, the rotational axis of a satellite coincides with the 
axis of the maximum moment of inertia of the satellite and is orthogonal to the orbital plane. 

Another qualitative kind of rotation known from observations is fast regular rotation. 
There are seven satellites that are known to rotate so (Maris et al., 2001, 2007; Grav et al., 
2003; Bauer et al., 2004; http://solarsystem.nasa.gov/planets/): Himalia (J6), Elara (J7), 
Phoebe (S9), Caliban (U16), Sycorax (U17), Prospero (U18), and Nereid (N2); all of them 
are irregular satellites. These satellites, apparently, are tidally unevolved. 

A third observationally discovered qualitative kind of rotation is chaotic tumbling. Wis- 
dom et al. (1984) and Wisdom (1987) demonstrated theoretically that a planetary satellite of 
irregular shape in an elliptic orbit could rotate in a chaotic, unpredictable way. They found 
that a unique (at that time) probable candidate for the chaotic rotation, due to a pronounced 
shape asymmetry and significant orbital eccentricity, was Hyperion (S7). Besides, it has a 
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small enough theoretical timescale of tidal deceleration of rotation from a primordial rota- 
tion state. Later on, a direct modelling of its observed light curves (Klavetter, 1989; Black 
et al., 1995; Devyatkin et al., 2002) confirmed the chaotic character of Hyperion's rotation. 
Recent direct imaging from the CASSINI spacecraft supports these conclusions (Thomas 
et al., 2007). 

It was found in a theoretical research (Kouprianov and Shevchenko, 2005) that two other 
Saturnian satellites, Prometheus (S16) and Pandora (S17), could also rotate chaotically (see 
also Melnikov and Shevchenko (2008)). Contrary to the case of Hyperion, possible chaos in 
rotation of these two satellites is due to fine-tuning of the dynamical and physical parameters 
rather than to a large extent of a chaotic zone in the rotational phase space. 

We see that the satellites spinning fast or tumbling chaotically are a definite minority 
among the satellites with known rotation states. However, the observed dominance of syn- 
chronous behaviour might be a selection effect, exaggerating the abundance of the mode 
typical for big satellites. This is most probable. Peale (1977) showed on the basis of tidal 
despinning timescale arguments that the majority of the irregular satellites are expected to 
reside close to their initial (fast) rotation states. 

A lot of new satellites has been discovered during last years. Now the total number of 
satellites exceeds 160 (see http://solarsystem.nasa.gov/planets/). The rotation states for 
the majority of them are not known. In what follows, we investigate the problem of typical 
rotation states among all known satellites. 

We consider the motion of a satellite with respect to its mass centre under the following 
assumptions. The satellite is a nonspherical rigid body moving in a fixed elliptic orbit about 
a planet. We consider the planet to be a fixed gravitating point. The shape of the satellite is 
described by a triaxial ellipsoid with the principal semiaxes a > b > c and the corresponding 
principal central moments of inertia A < B < C . The dynamics of the relative motion in the 
planar problem (i.e., when the satellite rotates/librates in the orbital plane) are determined 
by the two parameters: u = \J3(B — A)/C, characterizing the dynamical asymmetry of 
the satellite, and e, the eccentricity of its orbit. Under the given assumptions, the planar 
rotational/librational motion of a satellite in the gravitational field of the planet is described 
by the Beletsky equation (Beletsky, 1965): 

d 2 d9 
(1 + e cos /) — 2 ~ 2e sin / — + w 2 sin 9 cos 9 = 2e sin /, 
d/ df 

where / is the true anomaly, 9 is the angle between the axis of the minimum principal central 
moment of inertia of the satellite and the "planet - satellite" radius vector. 

As follows from an analysis of the Beletsky equation (see Melnikov and Shevchenko (2000) 
and references therein), for a satellite in an eccentric orbit, at definite values of the inertial 
parameters, synchronous resonance can have two centres in spin-orbit phase space; in other 
words, two different synchronous resonances, stable in the planar rotation problem, can 
exist. Consider a section, defined at the orbit pericentre, of the spin-orbit phase space. At 
ujq = 0, there exists a sole centre of synchronous resonance with coordinates 9 = mod it, 
d9/dt = 1. If the eccentricity is non-zero, upon increasing the value of u , the resonance 
centre moves down the d9/dt axis, and at a definite value of u (e. g., for e = 0.1 this value 
is ~ 1.26) another synchronous resonance appears. Following (Melnikov and Shevchenko, 
2000), we call the former synchronous resonance (emerging at zero value of ojq) the alpha 
mode, and the latter one — the beta mode of synchronous resonance. Upon increasing the 
(Jo parameter, the alpha and beta modes coexist over some limited interval of cuq (the extent 
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Figure 1: Location of the satellites with known radii in the u co -e" diagram. 



of this interval depends on the orbital eccentricity), and in the section there are two distinct 
resonance centres situated at one and the same value of the satellite's orientation angle. 
Such a phenomenon takes place for Amalthea (J5) (Melnikov and Shevchenko, 1998, 2000). 
On further increasing the Uq parameter, at some value of ojq the alpha resonance disappears, 
i.e., it becomes unstable in the planar problem, and only the beta resonance remains. 

The "wo-e" stability diagram is presented in Fig. 1. Theoretical boundaries of the zones 
of existence (i.e., stability in the planar problem) of synchronous resonances are drawn in 
accordance with (Melnikov, 2001). Regions marked by "la" and "lb" are the domains of sole 
existence of alpha resonance, "II" is the domain of sole existence of beta resonance, "III" is 
the domain of coexistence of alpha and beta resonances, "IV" is the domain of coexistence 
of alpha and period-doubling bifurcation modes of alpha resonance, "V" is the domain of 
non-existence of any 1:1 synchronous resonance, "VI" is the domain of sole existence of 
period-doubling bifurcation modes of alpha resonance. 

The solid circles in Fig. 1 represent the satellites with known uq. The open circles 
represent the satellites with the ojq parameter determined by means of an approximation of 
the observed dependence of uo on the satellite size r, accomplished following an approach by 
(Melnikov and Shevchenko, 2007). In total, the data on sizes and orbital eccentricities are 
available for 145 satellites (Karkoschka, 2003; Sheppard and Jewitt, 2003; Sheppard, Jewitt, 
and Kleyna, 2005, 2006; Porco et al, 2007; Thomas et al., 2007; http://ssd.jpl.nasa.gov/); 
so, there are 145 "observational points" in the stability diagram u UQ-e" in Fig. 1. The 
horizontal bars indicate three-sigma errors in estimating cu - They are all set to be equal to 
the limiting maximum value 0.21. 

From the constructed diagram we find that 73 objects are situated in domain V, and 
in domain lb 12 objects are situated higher than Hyperion (a sole solid circle in domain 
lb), while in domain lb there are 15 objects in total. Synchronous state of rotation does 
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not exist in domain V. For the majority of satellites in domain lb (namely, for those that 
are situated higher than Hyperion) synchronous rotation is highly probable to be attitude 
unstable as in the case of Hyperion. So, 73 satellites in domain V and 12 satellites in domain 
lb rotate either regularly and much faster than synchronously (those tidally unevolved) or 
chaotically (those tidally evolved). Summing up the objects, we see that a major part (at 
least 85 objects) of all satellites with unknown rotation states (132 objects), i.e., at least 
64%, cannot rotate synchronously. 

In summary, though the majority of planetary satellites with known rotation states ro- 
tates synchronously (facing one side towards the planet, like the Moon), a significant part 
(at least 64%) of all satellites with unknown rotation states cannot rotate synchronously. 
The reason is that no stable synchronous 1:1 spin-orbit state exists for these bodies, as our 
analysis of the satellites location on the "wo~e" stability diagram demonstrates. They rotate 
either regularly and much faster than synchronously (those tidally unevolved) or chaotically 
(tidally evolved objects or captured slow rotators). 

With the advent of new observational tools, more and more satellites are being discov- 
ered. Since they are all small, they are all irregularly shaped (Kouprianov and Shevchenko, 
2006). Besides, the newly discovered objects typically move in strongly eccentric orbits (see 
http://ssd.jpl.nasa.gov/; Sheppard and Jewitt (2003)). Therefore these new small satellites 
are expected to be located mostly in domain V of the u uo-e" stability diagram. Consequently, 
either fast regular rotation (most probable) or chaotic tumbling (much less probable), but 
not the ordinary synchronous 1:1 spin-orbit state, can be a typical rotation state for the 
newly discovered planetary satellites. 
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